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Abstract 

The behavior of polyampholytes near a charged planar surface is studied by means of Monte 
Carlo simulations. The investigated polyampholytes are overall electrically neutral and made up 
of oppositely charged units (called blocks) that are highly charged and of the the same length. 
The influence of block length and substrate's surface-charge-density on the adsorption behavior is 
addressed. A detailed structural study, including local monomer concentration, monomer mean 
height, transversal chain size, interface-bond orientation correlation, is provided. It is demonstrated 
that adsorption is favored for long enough blocks and/or high enough Coulomb interface-ion cou- 
plings. By explicitly measuring the chain size in the bulk, it is shown that the charged interface 
induces either a swelling or a shrinkage of the transversal dimension of the chain depending, in a 
non trivial manner, on the block length. 
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I. INTRODUCTION 



Polyelectrolytes (PEs, i.e., charged polymers) containing oppositely charged functional 
groups (blocks) are called polyampholytes (PAs). A well known example is provided by 
proteins (in particular gelatin). This class of material can be advantageously used as adsor- 
bate and like that offer many applications in the every day life, such as paper production, 
dewatering photographic films etc. 

In the last decade, many experimentalist researchers have investigated the adsorption 
behavior of PAs at charged substrates. For example, Neyret et al. JjJ found that PAs with 
a negative net charge can adsorb onto negatively charged latex particles in accordance with 
the theoretical prediction of Joanny More recently, using silicon substrates Mahltig et 
al. 5I 0] studied the influence of the pH on the PA adsorption properties, where the net 
charge of the PAs can be tuned in a well controlled manner. 

From a theoretical viewpoint, the behavior of PAs near a charged surface has been studied 
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by various authors |4|, H la l| • Joanny |4j showed that PAs with non-zero net charge can 
adsorb on a surface that carries the same sign of electric charge for two distinct regimes: 
high and low screening strengths. Later, Dobrynin, Rubinstein and Joanny jjj, using scaling 
law arguments, sketched typical adsorption regimes for salt-free environments. Based on a 



similar approach to that employed in Ref. [9|], scaling theories were developed to elucidate 
;he adsorption process of a single PA chain |7[ and many PA chains on a charged sphere 
s|. Shusharina and Linse ljj examined grafted diblock PAs at uncharged surfaces by a 
mean-field lattice theory and characterized the brush structure. 



As far as computer simulations are concerned, only a few studies 



12 1 were devoted 



to the problem of PA adsorption. Akinchina et al. [12j performed Monte Carlo (MC) 
simulations on diblock PAs at uncharged spherical particles. The brush structure as a 
function of the charge ratio of the two blocks was examined 3|. Closely related to our 
problem, MC simulations were carried out by McNamara et al. to characterize the 
molecular pattern recognition setting in upon adsorption of sequenced PEs (i.e., charge- 
ordered PAs or non-random PAs) to planar patterned surfaces. Using screened Coulomb 
interactions, they concluded that chain-entropy prevents efficient pattern recognition jll|. 
However, a detailed structural analysis and especially the monomer concentration profiles 
were not provided in Ref. [ll] so that a clear description of the adsorption process of non- 
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random PAs onto charged surfaces is still missing. The fact that only patterned surfaces 



(checkerboards and stripes) were considered in Ref. ll( makes also the process of adsorption 
more specific and complicated. 

The goal of this paper is to provide the mechanisms governing the adsorption of sequenced 
PEs (with a zero net charge) containing highly charged blocks onto charged planar surfaces. 
One can already anticipate and say that the overall chain adsorption properties are influenced 
by the following driving forces: 

(i) Chain-entropy tends to keep the chain in the solution. 

(ii) Electrostatic block-block correlations tend to collapse the chain into a globular form. 

(iii) A strong electrostatic interface-ion coupling favors chain adsorption. 

We limit our study to the most simple situation where a single PA chain, whose blocks are 
equisized and carry the same magnitude of electric charge, interacts with a charged planar 
interface. That way, it should be easier to identify the relevant mechanics setting in PA 
chain adsorption. Computer simulations are especially useful in this regime of charge since 
most analytical theories DDDO . are only suitable for weakly charged blocks. Our article 
is organized as follows: The simulation model is detailed in Sec. HU Results are presented 
in Sec. Hill and concluding remarks are provided in Sec. IIV1 



II. MODEL AND PARAMETERS 



The model system under consideration is similar to that recently investigated for the 



adsorption of fully positively charged PEs 



13, 



14j . Within the framework of the primitive 



model we consider a regular polyampholyte near a charged hard wall with an implicit solvent 
located at z > 0) of relative permittivity e r ~ 78. To avoid the appearance of image forces 
14j we suppose that the substrate (located at z < 0) below the interface (at z = 0) has the 
same dielectric constant. 

The negative bare surface-charge density of the substrate's interface is — 0oe, where 
e = 1.60219 x 10 _19 C is the (positive) elementary charge and o~q is the number of charges 
(— e) per unit area. The latter is always electrically compensated by its accompanying mono- 
valent counterions of charge Z + e (i.e., cations with Z + = +1) and diameter a. The block 
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PE corresponding to a perfectly regular polyampholyte is made up of 2iV& alternating posi- 
tively/negatively charged blocks: Each block contains M monovalent charged monomers of 
the same sign of charge with a (typical) diameter a. Hence by construction the PA chain, 
containing altogether N m = 2N^M charged monomers, is globally electroneutral. The cor- 
responding counterions [i. e., monovalent anions [Z- = —1) and cations] are also taken 
explicitly into account. Thereby, all the constitutive microions are monovalent and mono- 
sized with diameter a. All those particles are immersed in a rectangular L x L x r box. 
Periodic boundary conditions are applied in the (x, y) directions, whereas hard walls are 
present at z = (location of the charged interface) and z = r (location of an uncharged 
wall). It is to say that we work in the framework of the primitive cell model. 
The total energy of interaction of the system can be written as 



where the first (single) sum stems from the interaction between a microion % and the charged 
interface, and the second (double) sum stems from the pair interaction between microions % 
and j with = |r* — rj\. All these contributions to U to t in Eq. (JTJ) are described in detail 
below. 

Excluded volume interactions are modeled via a hardcore potential [3] defined as follows 




(1) 



+ i Uhs ( ri ^ + U coul{rij) + U FEN Eijij) + U LJ (rij)] , 



i,i<j 




0, for Tij > a 
oo, for < a 



(2) 



for the microion-microion one, and 




oo 



for a/2 < Zi < r — a/2 



otherwise 
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for the interface- micro ion one. 



The pair electrostatic interaction between two microions i and j reads 




(4) 
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where Ib = /5e 2 /(47reoe r ) is the Bjerrum length corresponding to the distance at which two 
charges e interact with 1/(3 = ksT . The electrostatic potential of interaction between a 
microion % and the (uniformly) charged interface reads 
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where A 



2n i B(7Q is ^e so-called Gouy-Chapmann length. The latter has two physical 
meanings: (i) It is the desorption length for an isolated ion of charge e (initially at contact 
with an oppositely charged wall) yielding an energy penalty of ksT; (ii) It also corresponds 
(in the Poisson-Boltzmann theory) to the distance from the wall at which surface charges are 
half-compensated by their counterions. In order to characterize the electrostatic interface-ion 
coupling we consider the dimensionless "Moreira-Netz" parameter S = y ^(|. A modified 
Lekner sum was utilized to compute the electrostatic interactions with periodicity in two 
directions. To link our simulation parameters to experimental systems we choose a = 4.25 
A leading to the Bjerrum length of water Is = 1.68a = 7.14 A at T = 298K. 

The PA chain connectivity is modeled by a standard finite extendible nonlinear elastic 
(FENE) potential for good solvent, which reads 



U F ENE(r) 
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Rl 



for r < Rq 



for r > Rn 



(6) 



with k = 27ksT/a 2 and Ro = 1.5a. The excluded volume interaction between chain 
monomers is taken into account via a shifted and truncated Lennard- Jones (LJ) potential 
given by 



U LJ (r) 



4f 
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r 



+ e, for r < 2 1/6 a 



(7) 



0, for r > 2 1 / 6 a 

where e = ksT. These parameter values lead to an equilibrium bond length I = 0.98a. 

The equilibrium properties of our model system were obtained by using standard canoni- 
cal MC simulations following the Metropolis scheme. In detail, single particle (translational) 
moves were applied to all the counterions (i.e., anions and cations) and the monomers with 
an acceptance ratio of 50%. For the sake of efficiency of the spatial sampling, we also apply 
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TABLE I: List of key parameters with some fixed values. 



Parameters 

T = 298K room temperature 

(JqL 2 Numbers of charge (— e) on the substrate 

Z± = ±1 microion valence 

a = 4.25 A microion diameter 

Ib = 1.68a = 7.14 A Bjerrum length 

L = 35a (x, y)-box length 

t = 105a z-box length 

M number of monomer(s) per positive/negative block 

N m = 64 total number of monomers 



time to time "global" translational displacements for the PA chain, where the whole chain 
is translated. The same acceptance ratio of 50% was applied there. 

All the simulation parameters are gathered in Table H] The classes of the simulated 
systems can be found in Table ILT1 For each class of systems characterized by a value of S, 
we have systematically varied the block length M from 2 to 32 with intermediate values of 
M = 4,8, 16, leading altogether to a set of 30 simulation runs. The monomer concentration 
is set to c m = jf^ ps 4.976 x 10 _4 a~ 3 leading to a PA volume fraction = 4^c m (a/2) 3 m 
2.605 x 10~ 4 (see also Table UJ). The total length of a simulation run is set to 3 — 4 x 10 6 MC 
steps per particle. Typically, the first 10 6 MC steps were discarded, and about 2 — 3 x 10 6 
MC steps were used to perform measurements. 



III. RESULTS AND DISCUSSION 
A. Preamble: Bulk behavior 



Before evoking the problem of a PA at a charged interface, it is instructive to first 
consider the more simple bulk situation. To our knowledge, the influence of the length 
(M) of the charged block on the PA chain-size has not been systematically investigated 
in the literature. At all events, the understanding of the behavior of PA chains near a 
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TABLE II: Classes of simulated systems characterized by their value of 3 = -f-. The number 
of counterions (cations and anions) ensuring the overall electroneutrality of the system is not 
indicated. The symbol indicates bulk systems. 

3 ao£_ 

oc*) 

0.9265 64 

1.853 128 

3.706 256 

7.412 512 



charged wall necessitates a knowledge of the limiting bulk system. To this end, we provide 
additional simulation data for the bulk situation (where H = 0). Thereby, the PA chain 
and their counterions are placed in a spherical cell of radius R ce ii = 31.31a such that the 
monomer concentration remains identical to that studied in the slab geometry. The 

' l7TH cell 

central (A^ m /2) th monomer along the chain is held fixed and standard single particle moves 
as described above are applied. The duration of the simulation is set to 10 8 MCS leading to 
excellent statistics. 

To characterize the chain conformation the radius of gyration defined as 

/ R ^\ J. yj(( rj - rom f) ( 8 ) 

Jv m , =1 

(where r cm designates the center of mass of the chain) was monitored. The results are 
presented in Fig. ^ Simulations show unambiguously that the chain shrinks with increasing 
block length M. This behavior is the result of an increasing cohesion between oppositely 
charged blocks with growing block charge. Similarly to what happens in an ionic crystal, 
upon increasing M one concomitantly enlarges the coordination number (i.e., the average 
number of neighbors of opposite charge). For longer chains, it is expected that the chain 
size should saturate from a certain threshold of M due to excluded volume effects. Some 
support to this latter statement seems to be provided by our simulations where the variation 
of the chain size is significantly slower at large values of M (see Fig. Q). Those preliminary 
findings will be useful to discuss our forthcoming results concerning the PA chain behavior 
near a charged interface. 
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FIG. 1: Bulk mean square radius of gyration (R^ 1 ^ ) as a function of the block length M. 
B. Behavior at weak interface-ion Coulomb coupling 

In this part, we focus on the situation where the interface-ion Coulomb coupling is weak 
(i.e., essentially 5 = 0.9265). Since the interface is (negatively) charged and the chain is 
flexible, different behaviors are expected for the positively and negatively charged blocks, 
in contrast to what occurs in the bulk. It is important to have in mind that the role of 
the flexibility of the chain is crucial here, since the stiff version of our PA would lead to a 
conformation where the rod-like chain lies parallel to the normal z-direction of the interface 
(with one end touching the interface) in a minimal energy configuration. 

To characterize the PA chain adsorption/depletion aspects, we have computed the 
monomer density profiles n±(z) defined as 

N m /2 

n ± (z)=J2 (*(*-** <±} ))' (9) 

i=l 

and normalized as follows 

m— a/2 nr 

/ n ± (z)L 2 dz = ^, (10) 

J a/2 1 

where +(— ) applied to positively (negatively) charged monomers. Our results are displayed 
in Fig. |2] and some microstructure snapshots are reported in Fig. El From a general 
viewpoint, the profiles of n+(z) and n_(z) in Fig. El are quite similar. The reason for 
that observation is (i) the connectivity between oppositely charged blocks and (ii) the weak 
Coulomb interface-block coupling. Those density profiles n±(z) exhibit a peak from M > 4 
which is a signature of PA adsorption (see Fig. |2J). The disappearance of a marked peak 
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FIG. 2: Density profiles n±(zf at S = 0.9265 for various values of M. 



in can also be used as a criterion 17[ for adsorption/depletion transition in our 

simulation. Such a criterion was utilized by Shafir el al. in their theoretical study of 
adsorption and depletion of polyelectrolytes from oppositely charged surfaces Near 
contact (z/a < 1.5) negatively charged monomers are strongly repelled from the wall due 




FIG. 3: Microstructure snapshots at H = 0.9265 for block lengths M = 1, 16,32 as indicated. The 
red (white) monomers correspond to positive (negative) charges. The little counterions are omitted 
for clarity. 
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to chain-entropy loss and wall-block electrostatic repulsion [see n_(z) in Fig. The 
positively charged monomers experience a short range repulsion as well but this time only 
due to chain-entropy effects [see n + (z) in Fig. This feature is well known for neutral 
chains Q|, and was also reported by several authors for charged ones 13, 18, 20, 21 1. The 
fact that the peak in n + (z) is considerably broadened (here at M > 4, see Fig. |2J) suggests 
that the adsorption is rather weak. This statement will be justified later where larger H 
coupling parameters are considered. It could be instructive to compare the contact value 
of n + (z — > a/2) to that obtain in the case the where only the substrate's counterions are 

n 

present 22] [i.e., the well known and exact contact value from the Gouy-Chapmann theory 
n^ c ^ = (2-kIbX 2 )^ 1 = Tjfp"]- Introducing the dimensionless quantity v defined as 

_ n + {z -» a/2) _ -> a/2) 

Fig. El shows that a typical value of v in the adsorption regime is roughly of the order of 
v ~ ° ° 2 ° 8 2 8 w 0.07. This low value confirms that the adsorption at H = 0.9265 is (very) 
weak. Close to the wall (z/a < 4), Fig. 121 indicates that n± increases systematically with 
M if the special diblock case (M = 32) is ignored. The non-trivial behavior at M = 32 is 
a result of a subtle competition between the interface-block interaction and the block-block 
correlations. The latter favor non-flat globular structures as well illustrated in Fig. |3] for 
M = 32. That effect can in turn lower the degree of adsorption and explain the non-trivial 
adsorption behavior reported in Fig. [21 at high N m . 

To further characterize the PA behavior at charged interfaces, we have measured the 
reduced mean height (h±) defined as 



r-o/2 



z 



(h±) = / n±{z)-L 2 dz, (12) 

Ja/2 a 

where +(— ) applied to positively (negatively) charged monomers. Within this definition 
the mean height of the center of mass, (h crn ), corresponds to (h cm ) = _ as 

a function of M (at S = 0.9265) are reported in Fig. HJ One can distinguish two typical 
regimes (see Fig. 0J: (i) At short charged blocks (M = 1,2) the values of (h±) are such that 



3- 



(h±) pa j- (within the statistical uncertainties) so that the chain is fully delocalized 
(ii) For longer blocks (M > 4), the mean heights decay by roughly one order of magnitude 
indicating an adsorption behavior that is consistent with the density profiles from Fig. El 
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FIG. 4: Reduced mean heights (h±) as a function of block length M 

We now would like to address the problem of desorption. The term "desorption" generally 
means that some adsorbed substance gets released, which implicitly involves a dynamic 
process. For our present study, this latter feature can be best captured by monitoring the 
height of the monomers as a function of time. For that issue, we define the reduced height 
of positive charges, h+(t), as follows 

= -H^T- < 13 > 

where the symbol "+" in the sum means that only positively charged monomers are counted. 
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FIG. 5: Reduced height h+(t) as a function of time for different parameters as indicated, to is a 
reference time. 
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The results for h + (t) for different set of parameters are shown in Fig. El The panels (a), 
(b) and (c) from Fig. El correspond to the lowest value of H = 0.9265 whereas panel (d) 
corresponds to a higher value of H = 1.853. The case of weakest interface-block coupling [i.e., 
H = 0.9265 and M — 1, see Fig. EJa)] shows large fluctuations of amplitude } l ^ nax '> _ } l ^ nm ') ^ 
r/a that are a dynamic signature of the derealization, in agreement with the value of the 
static quantity (h+) ~ j- found in Fig. 01 for M = 1. Upon increasing the interface-block 
coupling (i.e., increasing M and/or S), Fig. El indicates that the probability of escaping 
from the adsorbed state gets gradually weaker, as intuitively expected. The fact that at 
sufficiently high interface-block coupling the desorption processes become rare events, may in 
some cases considerably affect the quality of the statistical averaging as observed in Fig. HJfor 
M — 8. Consequently, a perfect equilibration is not always guaranteed for some parameters. 
Experimentalists are faced to the same problem when dealing with adsorption/desorption 
processes. This being said, we are confident that our simulations are reliable enough to 
capture the main qualitative features of the behavior of non-random PA chains at charged 
surfaces. 

C. Behavior at intermediate and large interface-ion Coulomb coupling 

In this part, larger interface-ion Coulomb couplings are considered corresponding to 
H = 1.853,3.706 and 7.412 (see Table ITT|) with varying block lengths M. Many observ- 
ables are going to be analyzed in detail and those can be classified into two main categories: 
(i) monomer concentration profiles and related quantities and (ii) chain conformation prop- 
erties. 

1. Monomer distribution 

The monomer density profiles are depicted in Fig. El and some corresponding microstruc- 
ture snapshots are sketched in Fig. Let us first discuss the behavior of the positively 
charged monomers which is relevant for the adsorption, since without those no adsorption 
at all would occur. At H = 1.853, the density profiles of the positively charged monomers 
exhibit similarity to that at H = 0.9265 (compare Fig. Eland Fig. |2J, where especially 
the repulsive chain-entropy effects are present. An adsorption behavior is detected here for 
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FIG. 6: Density profiles n±(z) 3 for different block lengths M. The values of 3 increase from top 
to bottom: (a) 3 = 1.853 (b) 3 = 3.706 (c) 3 = 7.412. 

M > 4 [see Fig. IH^a)]- This onset of adsorption can also be nicely visualized by compar- 
ing the snapshot of Fig. (at M = 4 and S = 1.853) to that of Fig. El (at M = 1 and 
H = 0.9265) corresponding to the non-adsorbing case. However, at given M, more posi- 
tively charged monomers can be found near the oppositely charged wall at H = 1.853 than 
at 5 = 0.9265 [roughly by a factor 2, compare Fig. Efa) and Fig. At higher wall-ion 
couplings (i.e., H = 3.706 and H = 7.412), the scenario gets qualitatively different (see Fig. 
If)}, where the n + (2;)-profiles are now monotonic in contrast to what was found at 5 = 0.9265 
(see Fig. EJ) and E = 1.853. Physically, this means that the electrostatic effects overcome 
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FIG. 7: Microstructure snapshots at different values of S and M as indicated. Note that S increases 
from top to bottom. The red (white) monomers correspond to positive (negative) charges. The 
little counterions are omitted for clarity. 
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the chain-entropy ones at sufficiently large S. This feature can also be seen as signature 
of strong adsorption. A closer look at Fig. E] reveals another interesting behavior, namely 
at our largest value of H = 7.412, the density n + (z) near contact increases systematically 
with growing block length M, in contrast to what happens at lower 5 (see also Fig. |2J). 
This phenomenon can be explained by the strong wall-block electrostatic interaction that 
overcomes the cohesive block-block correlations. Those both effects grow with M in a non 
trivial way. 

As far as the local concentration of negatively charged monomers n_ (z) is concerned, Fig. 
El indicates a complicated behavior as a function of S and M in the adsorption regime (here 
M > 4). More explicitly, at H = 1.853 the maxima reached by ri-(z) at long enough block 
length (here M > 8) are smaller than those of n + (z) [see Fig. E^a)], in contrast to what 
was reported at S = 0.9265 (see Fig. |2J). This behavior is due to the stronger wall-block 
correlation existing at H = 1.853. It is to say that the positively charged blocks get more 
attracted towards the negatively charged wall at higher S and concomitantly the negatively 
ones get more repelled from it. The density profile Ti-(z) at H = 3.706 is qualitatively 
similar to that obtained at H = 1.853 (see Fig. EI), and a truly different behavior is found at 
5 = 7.412. There, in the strong adsorption regime (i.e., H = 7.412 with M > 4), Fig. Et c ) 
interestingly indicates that only chains with intermediate block length M = 4, 8 can exhibit 
(significant) maxima in n_(z). The physical interpretation of this feature is that at high M 
the repulsion between the wall and the negatively charged blocks is so strong that, despite 
of the chain connectivity, the negatively charged get strongly repelled from the wall leading 
to formation of loops and tails. This mechanism explains also why n_(z), at H = 7.412, is 
higher for M = 4 than for M = 8 in Fig. E^c). All those relevant findings at H = 7.412 are 
convincingly illustrated in Fig. 0(at H = 7.412) where a "caterpillar" structure is developed 
for M = 8, whereas a "L" structure rises for M = 32. Note that the former structure was also 
theoretically reported by Dobrynin et al. [9] for symmetric PAs that was referred to as the 
"pancake regime" . On the other hand, the "L" structure deserves some further comments. 
Naively one could think that solely the strong electrostatic interface-block interactions are 



sufficient to explain this structure. Basic electrostatics [24 1 show that this idea is wrong, and 
one must take into account the energy gain upon approaching the two oppositely charged 
blocks when passing from a stretched configuration (perpendicular to the interface) to the 
"L" structure. Hence, the "L" structure is again the result of interface-block interactions 
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FIG. 8: Reduced mean heights {h±) as a function of block length M for various values of 5 
increasing from top to bottom: (a) E = 1.853 (b) E = 3.706 (c) E = 7.412. 

and block-block correlations. 

To gain further insight into the monomer distribution, we have plotted the mean heights 
(h±) as a function of M in Fig. |S| For all values of S, a delocalized state is reported at M = 1 
as signaled by the value of (h±) ~ r/2 (see Fig. |HJ). That non-adsorption at M = 1 is due to 
the very weak polarizability of the chain there. The problem of chain polarization [see Eq. 
(I14J) ] through the electric field stemming from the interface will be discussed in more detail 

— _ ° 2 

later. Given that our highest value S = 7.412 would correspond to |cr e| ~ le/(43A ) (see 
also Table HJfor the chosen values of system parameters), which represents already the upper 
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limit of experimentally accessible surface charge densities, we expect that no adsorption 
should occur in reasonable experimental conditions when M — 1. At M = 2, only the case 
H = 1.853 remains clearly delocalized (see Fig. [SJ). To rigorously assess the problem of 
localization/delocalization for the other remaining systems (especially for M = 2 at higher 
H), a more appropriate investigation involving a systematic variation of the box height r is 
required. Thereby if (h±) — > oo while r — > oo, then the chain is said delocalized, otherwise 
it is localized. Noticing that (h + ) — represents the (reduced) dipolar moment (p) of the 
chain given by 

0,> = /gz*\=jv m Mz<M, (14) 

a closer look at Fig. |H] reveals (at given M > 4) that the adsorbed chain gets more and 
more polarized with increasing S. The chain polarization is particularly vivid at high wall- 
block couplings where it increases with M [see Fig. Et c )]- Concomitantly, this M-enhanced 
polarization effect yields a strong non-monotonic behavior for (hJ) at S = 7.412 [see Fig. 
Etc)]. Nonetheless, we stress the point that, in general for intermediate and weak ion-wall 
couplings S, the polarization behavior as a function of M is non-trivial since (p) is also 
sensitive to the block-block correlations. 



2. Conformational properties 



A suitable quantity assessing the conformation of the chain is provided by the z- 
component {R 2 gz ) of the squared radius of gyration of the chain that is defined as follows 

W= ArE((*-«°») a >. (15) 

where z cm designates the z-component of the center of mass of the chain. The square root 



(Rg Z ) can also be seen as the thickness of the chain in the adsorbed state. It is instructive 
to compare those values to those obtained in the bulk case (see Fig. QJ. To this end, we 
introduce the following dimensionless quantity 

V = , (16) 

(R^ lk) ) 

where (R g b ^ lk ^ ) = ^(R g bulk ^ ). Thereby, rj < 1 signals an interface- induced compression 



whereas rj > 1 indicates an interface-induced stretching. Results for {R gz ) are presented in 
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FIG. 9: Mean square radii of gyration (Rg Z ) (slab) and (R^ lk ^ } (bulk) as a function of block 
length M for various values of 3 increasing from top to bottom: (a) 3 = 1.853 (b) 3 = 3.706 (c) 
3 = 7.412. 

Fig. El Despite of the rather large error bars appearing for the data points corresponding 
to the depletion/weak adsorption regimes, our simulations show that there is a systematic 
compression/stretching transition when passing from M — 1 to M — 2. This non-trivial 
phenomenon can be rationalized by considering the stiff limiting case. Indeed in the situation 
of a stiff rod-like PA, the minimal energy configuration (ignoring the counterions) would 
correspond to a chain "standing" on the wall with the positively charged end touching the 
interface, as previously discussed. Thereby, it is straightforward to show that the energy 
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gain increases with M (i.e., the dipolar moment of the chain). Hence at finite temperature, 
the rod-like will adopt an average orientation that gradually tends to be parallel to the in- 
direction when M is enlarged. Now, this idea can be qualitatively applied to flexible chains 
as long as the block-block correlations are not too strong, which is qualitatively the case for 
M = 1,2 (see Fig. [IJ. Our scenario is perfectly illustrated by the snapshot in Fig. at 
H = 3.706 and M = 2. This stretching at M = 2 implies also larger values of (R 2 gz ) than at 
M — 1 (see Fig. EJ). Note that this structure was also theoretically reported by Dobrynin 
et al. I| for symmetric PAs that was referred to as the "pole regime" . 

For longer blocks the situation is more complicated since there is a subtle competition 
between block-block correlations and block-wall ones. At intermediate block lengths (4 < 
M < 8) the chain is (again) in the compression regime (see. Fig. EJ), whereas by further 
increasing M the chain gets stretched. A possible exception to this scenario is the diblock 
case (M = 32) at H = 1.853 where 77 ~ 1 due to the strong block-block correlations and 
the rather weak interface- ion coupling there (see. Fig. EJ). The "L" structure adopted by a 
diblock polymer chain at S = 7.412 (see Fig. EJ) leads to a strong stretching stemming from 
the negatively charged block (see. Fig. EJ). 

The chain conformation can be further characterized by evaluating the bond-wall orien- 
tation correlation. To quantify this order parameter we monitor the following second order 
Legendre polynomial 

2 



' -^bond -i 

£5 



S±(z) 
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(17) 



where b 



(±) 



f (±) 



(with r 



(±) 



„(±) 



(±), 



x\ 'e x + y\ 'e v + z\ 



stands for the bond 



vector linking only positively (+) charged monomers or only negatively (— ) charged 

monomers, N bond = JVm ^~ 1 ^ is the total number of "positive" (or "negative") bonds, and 

=(±) 



mm z. 



' bond 



) Z i+l 



!. Hence, S±(z) defined for M > 2 reaches the values — |, 0, and +1 
for bonds that are perpendicular to the z-axis, randomly oriented, and parallel to the z-axis, 
respectively. Results for S±(z) are sketched in Fig. El Let us first analyze S+(z) where 
bonds link positively charged monomers. At contact {z < a), it can be seen from Fig. El 
that the positive bonds typically tend to lie the more parallel to the wall the larger the block 
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FIG. 10: Average bond- wall orientation correlation function S±(z) for different block lengths M. 
The values of 5 increase from top to bottom: (a) S = 1.853 (b) S = 3.706 (c) E = 7.412. 

lengths, as signaled by negatively stronger values of S + (z). Moreover, the wall-bond degree 
of parallelism is enhanced when H is increased at contact (see Fig. 110)1 . Sufficiently far away 
from the wall, bonds get typically randomly oriented as signaled by a near value reached 
by S + (z) (see Fig. ITUj) . Nonetheless, in the regime of strong adsorption (i.e., S = 7.412 and 
for M > 4), our simulation data show some deviation from that rule [see Fig. flDT c)]. This 
is merely due to the fact that at strong adsorption, positively charged monomers lying "far 
away" from the wall constitute rare events and therefore our results are not reliable for such 
high wall separations. 
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Concerning the bonds linking negatively charged monomers, an inspection of Fig. ITU1 
reveals that the S- (z)-profiles can be seen qualitatively as the mirror reflexion of the S + (z)- 
profiles about the axis S±(z) = 0. Still on a qualitative level, this mirror symmetry reflects 
simply the fact that positive bonds tend to get parallel to the wall whereas negative ones 
tend to stand perpendicular to the wall. This scenario is especially vivid at high S = 7.412 
as can be seen in Fig. HUT c) and in the snapshots from Fig. [7| taken at H = 7.412. In 
particular, for the diblock case (M = 32), the snapshot [see Fig. dc)] reveals that the tail 
is stretched leading to a strong and long-ranged wall-bond orientation correlation with a 
typical correlation length that should be proportional to Ma. 

IV. CONCLUDING REMARKS 

We have presented a systematic set of MC simulations to address the behavior of a 
single flexible charge-ordered polyampholyte chain near a charged planar interface. The 
chain length was fixed and the constitutive blocks were highly charged with two consecutive 
blocks being oppositely charged. The influence of block length M and interface-ion Coulomb 
coupling H on the overall chain conformation and adsorption properties were analyzed. 
Several relevant quantities were studied such as the local monomer concentration, the mean 
height of the centers of the positive and negative charges, the perpendicular component of 
the mean squared radius of gyration, the interface-bond orientation correlation function, 
and illustrative microstructural snapshots. Based on those results, our main findings can be 
summarized as follows: 

• No adsorption occurs at monomer-sized blocks due to the weak chain polarization 
there. 

• Chain adsorption can be obtained by increasing the interface-block coupling, i. e., 
increasing M and/or H. 

• By taking the bulk reference, compression-stretching transition sets in when 
passing from monomer-sized blocks to dimer-sized ones. The same behavior applies 
to the adsorption regime when passing from intermediate values of M to high values 
of M. 

• Formations of tails and loops manifest only for strong interface-block couplings. 
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These findings can be possibly useful to provide some hints to recent experiments [5( 
where the (lateral) hydrodynamic diameter of adsorbed diblock polyampholitic structures 
on silicon substrates and that in the solution were measured by light scattering techniques. 



y, the qualitative behavior of rj [Eq. 



is experimentally specified (see Fig. 



Consequent . 

n 

4 m Ref. [5J). Upon varying the pH Mahltig et al. |5j were able to gradually tune the net 
charge (from negative to positive passing through the isoelectric point) of the polyampholyte 
and that of the interface, and observed a similar compression-stretching scenario to that of 
our simulations. 
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